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摘要 本文 研究 合 随 机 风电 的 机 组 组 合 问题 , 基于 分 布 鲁 棒 优化 提出 了 考虑 风 功 率 预 测 的 条 件 误差 
和 误差 时 空 相关 性 的 建 模 及 优化 方法 . 首先 , 为 了 量化 计算 风 功 率 的 预测 误差 , 用 Lasso (最 小 绝对 
值 收 伍 和 选择 算 子 ) 训练 得 到 和 鲁 棒 的 条 件 误 差 估 计 器 ; 为 了 获取 风 功 率 预 测 误差 的 时 空 关 联 信息 , 通 
过 无 偏 佑 计 从 历史 数据 得 到 误差 的 协 方差 矩阵 .用 所 得 的 条 件 误 差 和 协 方 差 矩 阵 构 造 了 改进 的 风 功 
率 预 测 误差 的 概率 分 布 模糊 集 ， 其 次 ,在 多 面体 支撑 集 上 建立 了 风 功 率 预 测 误 差 的 随机 量 , 构建 了 两 
阶段 分 布 鲁 棒 机 组 组 合 模型 ,并 提出 了 与 该 模型 等 效 的 混合 整数 半 正 定 规划 模型 . 再 次 , 提出 了 从 第 
二 阶段 半 正 定 优化 问题 辨识 极限 分 布 的 方法 ， 并 提出 了 利用 极限 分 布 求解 两 阶段 分 布 鲁 棒 机 组 组 合 
问题 的 高 效 割 平面 算法 . 最 后 , 进行 数值 实验 , 验证 了 所 提 模 型 在 应 对 风 功 率 预 测 误差 的 时 空 关 联 性 
方面 的 能 力 和 优势 , 验证 了 模型 决策 方案 的 经 济 性 和 重 棒 性 . 

Ki MAKER, 分 布 鲁 棒 优 化 , 极限 分 布 , 时 空 关联 性 , 风 功 率 预 测 


1 引言 


为 应 对 化 石 能 源 枯竭 问题 和 环境 污染 问题 ， 风 电 等 清洁 能 源 发 电 在 近 二 十 年 来 得 到 了 快速 应 用 
和 发 展 . 风电 在 助力 全 球 可 持续 发 展 的 同时 ,也 给 电力 系统 的 经 济 和 可 靠 运行 ,如 日 前 机 组 组 合 , 带 
来 了 新 的 挑战 [Qin et al.(2016),Lowery and O’Malley(2012)}. 这 一 挑战 的 根源 是 风 功 率 的 预测 精度 
BUR, 尤其 是 对 于 在 地 理 上 广泛 分 布 的 多 个 风电 场 . 为 了 高 效 地 消 纳 风电 并 规避 风电 的 负面 作用 , 机 
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组 组 合 问题 需要 精细 化 考虑 风电 的 随机 性 ， 以 更 好 的 开 停机 计划 来 保证 风 功 率 随 机 扰动 下 的 系统 功 
率 平衡 和 经 济 运行 . 

随机 规划 是 决策 问题 中 随机 变量 建 模 的 经 典 方法 ， 但 其 在 电力 系统 机 组 组 合 问题 的 应 用 中 显露 
出 多 方面 的 不 足 . 一 方面 , 风 功 率 预 测 误差 的 概率 分 布 实际 上 是 无 法 准确 估计 的 , 因此 随机 规划 是 建 
立 在 一 个 有 偏差 的 概率 分 布 之 上 ， 其 得 到 的 决策 方案 自然 是 有 一 定 误 导 性 的 [Delage and Ye(2010)]. 
另 一 方面 ， 随 机 规划 中 的 数学 期 望 无 法 准确 计算 ， 需 要 借助 蒙特 卡 洛 模拟 等 技术 来 估算 (或 者 称 样 
本 平均 估计 ) [Dupacova et al.(2003)], ， 而 场景 采样 过 程 必然 损失 精度 . 所 以 ， 随 机 规划 得 到 的 开 停 
机 方案 可 能 是 不 可 靠 的 ， 再 调度 中 可 能 需要 付出 昂贵 的 代价 . 鲁 棒 优化 是 风 功 率 不 确定 性 建 模 的 为 
一 种 方法 [Bertsimas et al.(2011), Zeng and Zhao(2013), Samani and Aminifar(2019)]. 通过 考虑 不 
确定 集 内 最 恶劣 的 场景 ， 鲁 棒 优 化 可 以 得 到 可 靠 性 极 高 的 决策 方案 . 基本 上 ， 重 棒 优 化 不 考虑 不 确 
定量 的 概率 分 布 和 详细 关联 性 信息 [Zhao and Zeng(2012), Bertsimas et al.(2013)], ， 而 这 些 正 是 风电 
的 关键 信息 . 近年 来 ， 文 献 [Lorca and Sun(2016)] 提出 了 计 及 可 再 生 能 源 发 电 时 空 关联 性 的 动态 不 
确定 集 ， 文 献 [Roldan et al.(2019)] 基于 协 方 差 矩 阵 提 出 了 计 及 不 确定 量 相 关 性 的 椭 球 型 不 确定 集 . 
然而 ， 这 些 文献 中 的 不 确定 集 并 不 包含 概率 分 布 信息 . 

精细 和 严谨 地 考虑 风电 的 随机 性 对 于 新 型 电力 系统 的 运行 是 非常 必要 的 .基于 上 面 的 讨论 和 
分 析 ， 本 文采 用 随机 规划 和 重 棒 优 化 之 外 的 另 一 种 优化 方法 来 对 风电 随机 性 进行 建 模 ， 从 而 提升 
机 组 组 合 方案 的 经 济 性 和 可 靠 性 ， 分 布 鲁 棒 优 化 在 电气 工程 学 科 是 一 种 非常 新 的 范式 [Zhao and 
Guan(2015), Wei et al.(2016), Duan et al.(2017), Xiong et al.(2017),Chen et al.(2018),Zhu et al.(2019), 
Babaei et al.(2019), Zheng et al.(2020)], 其 建 模 思路 是 找到 给 定 模糊 集中 最 恶劣 的 概率 分 布 , 并 依据 
该 分 布 制定 决策 [Delage and Ye(2010)]. 相 比 随机 规划 ， 分 布 鲁 棒 优 化 借助 随机 量 的 矩 信 息 [Delage 
and Ye(2010), Bertsimas et al.(2010)] 或 者 一 个 训练 数据 集 [Esfahani and Kuhn(2018), Zheng and 
Chen(2020)] 来 定义 概率 分 布 的 集合 ， 考虑 最 极端 概率 分 布 得 到 的 决策 相 比 随机 规划 的 决策 具有 更 
高 的 可 靠 性 . 而 相 比 鲁 棒 优 化 ， 分 布 鲁 棒 优 化 考虑 了 更 多 的 概率 和 统计 信息 ， 即 它 的 决策 过 程 量化 
地 考虑 所 有 随机 事件 ， 而 不 是 聚焦 在 极端 场景 上 所以， 分布 鲁 棒 优 化 得 到 的 决策 在 期 望 的 意义 上 
通常 具有 更 好 的 经 济 性 . 

当前 ， 应 用 分 布 鲁 棱 优 化 处 理 机 组 组 合 问题 的 文献 大 多 采用 基于 概率 密度 的 模糊 集 [Zhao and 
Guan(2015), Duan et al.(2017)] 和 基于 概率 分 布 距离 的 模糊 集 [Chen et al.(2018),Zhu et al.(2019)， 
Zheng and Chen(2020)]. 也 有 文献 采用 基于 和 矩 信 息 的 模糊 集 ， 但 基本 上 只 考虑 了 一 阶 矩 信息 [Xiong 
et al.(2017)]， 所 以 没有 办 法 对 风电 场 之 间 的 时 空 关联 性 进行 建 模 . 尽管 有 文献 考虑 了 部 分 的 二 阶 矩 
信息 ( 即 方差 ) [Babaei et al.(2019)]， 但 方差 并 不 包含 随机 量 之 间 的 协 变性 . 实际 上 ， 协 方差 反映 
了 多 个 风电 场 出 力 的 时 空 关 联 性 ， 其 携带 的 信息 有 助 于 排除 风 功 率 不 符合 实际 的 快速 波动 和 有 疏 坡 事 
件 [Qin et al.(2016), Cordova et al.(2018)] , 而 这 些 极端 事件 能 够 显著 影响 机 组 的 开 停机 计划 [Lowery 
and O’Malley(2012), Che et al.(2019)]. 

KF FE OMI 75 — 7 ae EA URE. A, Serbs ARO AEN] 
再 生 能 源 发 电 概 率 预测 方面 广泛 应 用 [Cordova et al.(2018), Bessa et al.(2012), Agoua et al.(2019)]. 
例如 ， 文 献 [Bessa et al.(2012)] 提出 了 条 件 性 概率 密度 估计 器 ， 实 现 具有 时 间 适 应 性 的 风 功率 预测 ; 
文献 [Agoua et al.(2019)] 采用 Lasso 来 实现 时 空 预测 模型 的 输入 变量 筛选 . 不 过 ， 目 前 已 有 的 预测 
方法 通常 忽略 了 预测 值 和 预测 误差 的 关联 性 . 实际 上 ， 由 气象 或 物理 方法 给 出 一 个 日 前 预测 值 之 后 ， 
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该 预测 值 的 误差 期 望 是 有 一 定 规律 ， 并 且 可 以 进行 统计 推理 的 . 所 以 ， 可 以 对 条 件 预 测 误差 进行 量 
化 分 析 ， 然 后 用 于 校正 日 前 预测 ， 从 而 构建 更 为 精确 的 模糊 集 . 
除了 建 模 之 外 ,两 阶段 分 布 鲁 棒 机 组 组 合 模型 在 求解 方面 也 极 具 挑 战 性 , 尤其 是 考虑 了 二 阶 矩 约 
束 之 后 . 已 有 的 求解 方法 利用 线性 决策 规则 [Xiong et al.(2017), Duan et al.(2017), Zhu et al.(2019)， 
Babaei et al.(2019)], Benders 分 解法 [Zhao and Guan(2015), Chen et al.(2018)] ， 或 者 延迟 约束 生 
成 法 [Wei et al.(2016), Chen et al.(2016), Zheng et al.(2020)]. 线性 决策 规则 作为 一 种 近似 方法 ， 会 
导致 决策 的 次 优 性 . Benders 分 解 方法 和 延迟 约束 生成 法 则 存在 计算 效率 低 的 问题 ， 可 能 无 法 有 效 
处 理 模 糊 集 含有 协 方差 矩阵 的 分 布 鲁 棒 机 组 组 合 问题 [Zheng et al.(2020)]. 

在 这 一 项 研究 中 ， 我 们 使 用 精细 的 模糊 集 构建 了 基于 和 矩 模糊 集 的 分 布 鲁 棒 机 组 组 合 模 型 ， 从 而 
对 风 功 率 预测 误差 的 相关 性 和 条 件 误 差 进 行 了 有 效 建 模 . 此 外 ， 提 出 了 高 效 的 计算 方法 . 本文 的 贡 


ee 献 可 以 总 结 为 以 下 两 点 
K 1 利用 Lasso 从 训练 数据 集 闪 识 风 功 率 的 条 件 误差 与 日 前 预测 值 的 关系 . 使 用 风 切 率 预测 误差 的 


条 件 期 望 和 协 方差 构建 模糊 集 ， 从 而 建立 了 含 改 进 模糊 集 的 分 布 鲁 棒 机 组 组 合 模型 ， 


2. 利用 锥 线性 规划 的 对 偶 理 论 、 单 阶段 鲁 棒 优 化 的 重 构 方法 、Schur 补 等 ， 将 含 多 面体 支撑 集 的 
分 布 鲁 棒 机 组 组 合 模型 的 第 二 阶段 问题 精确 地 转化 为 半 正 定 规划 问题 . 证 明 半 正定 规划 问题 求 
解 之 后 可 以 得 到 一 个 风 功 率 的 极限 分 布 . 基于 极限 分 布 ,提出 了 基于 原 空间 制 平面 的 算法 ， 该 
算法 相 比 已 有 算法 能 够 更 加 高 效 的 求解 两 阶段 分 布 鲁 棒 机 组 组 合 问题 . 


2 模型 构造 
i 21 ”分布 鲁 棒 机 组 组 合 模型 
.全 首先 给 出 本 文采 用 的 机 组 组 合 模型 . 该 模型 考虑 了 可 中 晰 负荷 、 随 机 风电 、 弃 风 等 要 素 ， 模 型 
全 具体 如 下 所 示 . 
min 》 X zosNLy + ug tSUg + vosSDy + Co(post) 
teT gEG 
和 
tET dEDI tET wEW 
+>》 (S (St + 865) 二》 sD,) (1) 
tET gEG deD 
S.t. vot — Tot 1—= Ugt— Vgt VO SG LET (2) 
t 
Ugr S Zg VO CUTE T (3) 


T=max{1,t—MU,+1} 
t 
Ugr tee YJEG,tET (4) 
T=max{1,t-—MD,+1} 


Pee, S pge S Pag YJEG,tET (5) 


chinaXiv:202208.00027v2 


— R} S Pot — Pee S R YgEG,t,t-1ET 
0< Sie < DR vdeD',teT 
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-Fi < X KSilpoetsSt — s9; )+ 


gEG 
2 2Dat—sas) — > nardat 
deD deDI 
ynet ELET 
wEW 
= 
gEG wEW 


>》 (Dat—shs)— >》 Oe WET. 


2 去 dep! 
T: 调度 时 段 的 集合 

9: 发 电机 组 的 集合 . 

C: 输电 线路 的 集合 . 

D: 所 有 负 蓓 节点 的 集合 . 

D1: D 的 子 集 ， 代 表 可 中 断 负荷 的 集合 
W: 风电 场 的 集合 

NL,: 发 电机 组 9 的 空 载 成 本 

SU,: 机 组 9 的 启动 成 本 

SDs: 机 组 9 的 停机 成 本 . 

CoC): 机 组 g 的 可 变 发 电 成 本 . 

OLC): 负荷 d 在 时 段 t 的 负荷 中 断 成 本 
CSO): 风电 场 w 的 弃 风 成 本 . 
CP?en(.): 松弛 变量 的 惩罚 费用 . 

MU,: 机 组 g 的 最 小 开机 时 间 . 

MD,: 机 组 g 的 最 小 停机 时 间 . 

Prins 机 组 9 的 最 小 技术 出 力 . 

Prax: 机 组 9 的 最 大 发 电能 
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© RE: 机 组 g H E/F NEIRE 


Dp: 负荷 d 在 时 段 二 可 中 断 量 . 


Fy: 输电 线路 ! 的 传输 容量 . 


Koy? 机 组 g 到 线路 ! 的 功率 转移 因子 . 


。 KD: 负荷 d 到 线路 1 的 功率 转移 因子 
。 KW: 风电 场 w 到 线路 ! 的 功率 转移 因子 . 


© Das: 负荷 d 在 时 段 t 的 需求 水 平 . 
。 Wor: 风电 场 w 在 时 段 t 的 日 前 发 电功率 预测 值 . 


。 zot: 指示 机 组 开机 的 0-1 状态 变量 . 


。 Uge: 指示 机 组 启动 的 0-1 状态 变量 . 


。 Uge: 指示 机 组 关 停 的 0-1 状态 变量 . 


。 pot 机 组 g ÆRE t 的 出 力 优化 值 . 


ðk: 负荷 d 在 时 段 寺 的 负荷 中 断 量 . 


Ov: 风电 场 w 在 时 段 t WFE. 


。 sE: 机 组 g 在 时 段 t 松弛 变量 的 正 / 负 分 量 . 
sR: 负荷 d TERTE: t 的 非 负 松弛 变量 
e Eoi 风电 场 w 在 时 段 t 的 预测 误差 ， 随 机 量 


目标 函数 (1) 考 虑 了 开 停 机 成 本 、 发 电 成 本 、 负 荷 中 断 和 弃 风 费用 ,以 及 松弛 变量 的 惩罚 费用 等 
各 项 费用 均 采 用 线性 的 成 本 函数 . 式 (2)-(4 为 机 组 的 状态 转移 方程 和 最 小 开 停 机 时 间 约 束 . 式 (5) 为 
机 组 的 发 电能 力 约束 ， 该 约束 耦合 了 机 组 启 停 0-1 变量 和 出 力 连 续 变量 ， 式 (6) 为 机 组 的 息 坡 约束 . 
式 (7)-(3) 限 定 了 负荷 中 断 量 和 弃 风 量 . 式 (9) 为 输电 线路 的 潮流 约束 . 式 (10) 为 系统 功率 平衡 约束 . 注 
意 约束 (9)-(10) 包 含 松 弛 变量 ， 用 以 保证 问题 的 可 行 性 . 具体 而 言 ， 每 一 台 发 电机 组 的 每 个 时 段 均 配 
置 了 自由 变量 ， 表 征 正 的 功率 注入 ST 或 者 负 的 功率 注入 s97 ， 而 每 个 负荷 阶段 均 配 置 了 切 负 荷 变 
E s9,. 松弛 变量 取 非 零 值 意 味 系统 在 正常 情况 下 无 法 维持 线路 潮流 不 越 限 ， 或 者 不 能 维持 系统 功 
率 平衡 ， 必 须 采 取 紧 急 措施 才能 保证 安全 运行 . 为 了 表示 非 正常 运行 情况 下 的 风险 代价 ， 式 (1) 的 最 
后 一 项 对 松弛 变量 进行 惩罚 ， 该 惩罚 系数 取 值 较 大 

在 问题 (1D)-(10) 中 ， 机 组 启 停 变 量 是 “此 时 此 地 ”类 型 的 决策 变量 ， 后 续 将 其 归 类 为 第 一 阶段 
变量 z， 连 续 型 的 调度 变量 可 以 根据 日 内 观察 到 的 实际 情况 进行 调整 ， 因 此 归 类 为 第 二 阶段 变量 
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y [Bertsimas et al.(2013)]. 随机 变量 记 为 向 量 €. 基于 变量 的 分 组 ， 可 以 得 到 问题 (1)-(10) 的 紧凑 形 
式 ， 具体 如 下 所 示 . 


min clz 二 coy (lla) 

TY 
s.t. Fe >f (11b) 
Y, (110) 
Az + By 2d (11d) 
(11e) 


lle 


HP, cr 和 cc 是 成 本 系数 . 式 (11b) 包 含 hs )-(4). 式 (11c) 是 第 二 阶段 变量 y 本 身 的 相关 约束 ， 包 
含 (6)-(7). 式 (11d) 是 耦合 约束 ,对 应 (5). 问题 (1)-(10) 中 既 包 含 第 二 阶段 连续 变量 又 包含 随机 量 的 约 
束 由 式 (116) 表 示 . 

问题 (1)-(10) 对 应 的 两 阶段 分 布 鲁 棒 优化 模型 考虑 模糊 集中 最 极端 的 概率 分 布 进行 优化 决策 , 其 
具体 形式 为 


min max c] æ + Er[Q(z,é)], (12) 


wEex P 


HP X = {a|Fa > f} 是 第 一 阶段 变量 的 可 行 域 . 式 (12) 的 最 优 值 函数 Q(x,é&) 是 根据 第 二 阶段 问 
题 定义 的 ， 即 
Q(x, £) = min cy, (13) 


其 中 Y = {yl(11c) 一 (11e)} 是 第 二 阶段 变量 的 可 行 域 . 模糊 集 P 的 定义 如 下 ， 


Saye Ee} (15) 


式 (14) 的 第 二 行 表示 风 功 率 预 测 误差 的 期 望 ee ee ate Maca 
的 交集 之 内 ， 即 0 Elé- EE-E 直观 地 讲 ， 风 功率 预测 误差 的 方差 和 协 方差 不 能 
经 验 数值 给 定 的 程度 更 “严重 ”， 


2.2 ”基于 数据 的 矩 估 计 

随机 向 量 & 的 一 阶 和 二 阶 矩 对 于 模糊 集 (14) 的 构建 十 分 关键 . 在 这 一 小 节 ， 我 们 将 介绍 & 和 允 
的 佑 计 方法 . 

通常 ， 风 功率 日 前 预测 采用 气象 或 物理 方法 ， 由 数值 天 气 预报 和 风机 功率 曲线 得 到 风 功 率 预 测 
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值 [Shokrzadeh et al.(2014)]. 在 执行 机 组 组 合计 算 之 前 ， 系 统 运行 机 构 已 经 掌握 了 风 功 率 的 预测 值 . 
举 个 例子 , 在 一 个 具有 2 个 风电 场 的 系统 中 ， 如 果 采 用 24 小 时 的 计划 模型 ， 则 在 机 组 组 合 建 模 时 已 
经 具备 了 维 数 为 N = 48 的 日 前 预测 向 量 w. 给 定 日 前 预测 ,预测 误差 的 条 件 期 望 通常 是 非 零 值 ， 其 
至 可 能 很 大 [Bessa et al.(2012)]. 我 们 将 与 向 量 © 的 第 ”个 元 素 相对 于 的 预测 误差 条 件 期 望 标记 为 
[En [co]. 这 里 所 描述 的 条 件 期 望 便 是 模糊 集 (14) 所 用 到 的 期 望 ， 即 Elem] = &. 为 了 推断 出 Ele, |] 
的 数值 ， 假 定 我 们 具备 一 组 历史 的 日 前 预测 和 实际 误差 的 数据 ， 这 组 数据 的 形式 为 


Dp, E { (GiGi); (Qi Ens), e. (5, En,s)} 


E (Es, s) 具备 如 下 线性 结构 


N 
Gis = Bn,o “Ir 人 + En,s s= 1, Aa S. (16) 
J 


式 (16) 中 ，3 表示 数据 集 的 规模 ，snws。 为 噪声 ，{Bw0,…, Bn,N} 是 需要 辨识 的 权重 . 上 述 线性 模型 可 
以 写成 矩阵 的 形式 ， 


En = QBn + En, (17) 
其 中 Ên E ends Ens » En 二 [En,1) iia el Bn = [Bn,o; vey Bn â = [1's a]; said): 下 一 
步 的 任务 是 找到 一 个 能 够 最 好 程度 拟 合 历史 数据 的 向 量 成 ， 基 于 该 向 量 我 们 可 以 从 新 的 样本 w 计 
算出 &，。( 即 给 定 一 个 新 的 日 前 预测 值 ， 可 以 计算 出 相应 的 条 件 误差 ). 


| 


式 (17) 中 的 向 量 Ge 可 以 经 由 回归 分 析 得 到 . 在 多 种 回归 分 析 方 法 中 ，Lasso 相 比 经 典 的 线性 回 
归 和 怜 回 归 有 一 定 的 优势 ， 因 此 本 文采 用 Lasso 来 求解 线性 系统 (17). Lasso 是 基于 对 拟 合 偏差 平方 
和 与 Ly 范 数 惩罚 的 最 小 化 实现 的 ， 即 


"E. 
a — a sllén oan OB,,||3 + Allalla- (18) 


Ly 范 数 使 得 Lasso 具备 了 自动 特征 筛选 的 功能 ， 可 以 防止 过 拟 合 ， 实 现 更 强 的 鲁 棒 性 . 参数 A 是正 
则 化 的 权重 稀疏 ， 该 参数 可 以 由 交叉 检验 来 选取 [Lv et al.(2013)]. 本 文 的 Lasso 将 由 经 典 的 坐标 下 
降 算 法 求解 [Friedman et al.(2007)]. 


根据 式 (17)-(18), 所 有 输入 变量 统一 用 于 估计 任意 一 个 输出 . 所 以 , 估计 器 计 及 了 相 邻 时 段 和 相 
邻 风 电场 的 特征 . 值得 注意 的 是 ， 在 本 文 的 框架 内 需要 求解 N 个 Lasso 模型 ， 通过 这 种 方式 反映 每 
个 & 与 % 之 间 的 不 同 映射 关系 [Bessa et al.(2012)]. 最 终 ， 基 于 训练 完成 的 N 个 估计 器 ， 一 阶 矩 
可 以 在 日 前 风 功 率 预测 值 的 基础 上 计算 得 到 : 


& = |e be n=, (19) 
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随机 向 量 的 二 阶 和 矩 〈 协 方差 矩阵 ) 由 无 偏 估计 器 得 到 [Delage and Ye(2010)]: 


a eo > 人 -过 -二 > E) T. (20) 


3 求解 方法 


3.1 求解 第 二 阶段 问题 


分 布 鲁 棒 机 组 组 合 模型 (12) 的 第 二 阶段 最 大 化 问题 是 一 个 无 穷 维 锥 线性 问题 ， 其 详细 模型 为 


Sn me | Weedede 

s.t. fe>0 VEEE 
fiat =1 :ho (21) 
J std =€ ch 


j: EET fede 2 S+ EE H, 


其 中 fe 是 决策 变量 , 代表 事件 E 处 的 概率 密度 .fe 是 一 个 半 正 定 矩 阵 , h 是 自由 向 量 , ho 是 自由 变 
量 . 这 三 个 变量 是 问题 的 约束 的 对 偶 变 量 . 问题 21 的 对 偶 问 题 为 [Bertsimas and Sethuraman(2000), 
Ereund(2004)] ， 


Z(ax) = i, tr(H'@)+h'E+ho (22a) 
s.t. €'HE+h'E+ho > Ql, E) VEEE, (22b) 


其 中 © = N+ €€", tr(X) RRM X 的 迹 ， 值得 注意 的 是 ， 问 题 (22) 和 问题 (21) 的 强 对 偶 性 的 
充分 条 件 是 这 两 个 问题 都 是 可 行 的 ， 且 © > 0 [Bertsimas and Sethuraman(2000), Wei et al.(2016)], 
而 这 两 个 条 件 是 成 立 的 . 一 方面 ， 原 问题 (21) 总 是 可 行 的 ， 因 为 其 目标 函数 中 的 值 函 数 Q(x,é&) 总 是 
可 行 和 有 界 的 (由 第 2.1 节 机 组 组 合 问题 的 定义 可 知 ); 另外 , 矩 问题 至 少 有 一 个 可 行 解 ， 该 解 是 在 期 
望 & 处 取 概 率 1 的 分 布 [Delage and Ye(2010)]. 另 一 方面 ， 协 方差 矩阵 X 在 实际 中 通常 不 会 是 0, 
因为 风 功 率 的 日 前 预测 不 可 能 在 整个 24 小 时 范围 内 都 是 准确 的 . 

给 约束 (1lc)-(1le) 赋 予 对 偶 变 量 A, y, v 之 后 ， 可 以 写 出 Ql, é) 的 对 偶 问题 ， 


Q(z, £) = max g'A (Ax — d)' w—(VE-Ww)'v 


s.t.ecc —(G'A+ B'p+U'v)=0,A>0,n>0, (23) 


《 
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基于 A, u, v 的 对 偶 问 题 ， 问 题 (22) 可 以 写成 [Wei et al.(2016)]， 


Z(a)= min, tr (H'O) +h'é+ ho (24a) 
s.t. €'HE+h'E+ho > Q (#,€)—-(VE)'H VEEB ict, (24b) 


HP Q7 (æ, €) = g' Ai- (Ax 一 qd) p+Ww'v; 是 问题 (23) 的 目标 函数 中 独立 于 的 部 分 ,{ 和 Xi, His vi} 


问题 (24) 仍 然 是 一 个 半 无 穷 规 划 问 题 ， 该 问题 包含 无 穷 多 个 以 & 为 指标 的 约束 ， 因 此 无 法 直接 
求解 . 在 文献 [Wei et al.(2016), Zheng et al.(2020)] 中 ， 基 于 三 = RIS 的 假设 ,约束 (24b) 可 以 用 一 
个 半 正 定 锥 等 效 表 征 . 然而 ， 支 撑 集 无 穷 大 的 假设 要 求 系统 能 够 应 对 大 量 非 现实 的 风 功 率 场景 ， 进 
而 可 能 导致 最 终 的 机 组 组 合 方案 不 合理 . 在 文献 [Chen et al.(2016), Delage and Ye(2010)] 中 ， 基 于 
支撑 集 为 椭 球 的 假设 ,约束 (24b) 可 以 用 S 引 理 等 效 转化 成 一 个 半 正 定 约束 . 然而 ， 该 方法 不 能 直接 
应 用 于 约束 (24b)， 因 为 这 里 采用 的 是 多 面体 支撑 集 . 为 此 ， 本 文 提出 一 种 新 的 方法 来 处 理 半 无 穷 维 
约束 (24b). 该 方法 的 出 发 点 是 将 问题 (24) 看 成 一 个 以 & 为 不 确定 量 的 单 阶段 鲁 棒 优化 问题 . 该 方法 
具体 如 下 所 示 . 

尽管 H 并 不 一 定 是 整定 矩阵 ， 我 们 首先 假设 A ~ 0. 这 一 假设 意味 着 (25) 中 的 最 小 化 问题 是 
一 个 含 仿 射 约束 的 凸 二 次 规划 . 所 以 ， 根据 强 对 侦 性 ， 约 束 (26) 的 第 i 行 (后 续 记 为 (25)) 可 以 等 效 
地 写成 下 式 ， 


min §' HE +(h! +v; V)E + ho — Q; (x, €) > 0 (25) 
dualizing 
ma. —1/4h;(u;)' H~*h,(u;) T hoi (wi) 之 0. (26) 


式 (26) 中 ，wzi, wii 是 《的 边界 约束 对 应 的 非 负 对 偶 变 量 ， 具 体 参 照 式 (15).， 向 量 hilu) 表示 h+ 
TV + Ug; 一 Ui, bre hoi (u;) 表示 ho 一 Qr (x, £) F uzt 十 “Le. 注意 到 式 (26) 中 的 max 算 
子 可 以 移 除 ， 式 (26) 可 以 通过 Schur 补 进 一 步 写 成 


| H  1/2hj(us) 


| = 0, wii, Uzi 2 0. (27) 


对 于 五 =0 的 情况 , 式 (26) 将 退化 成 线性 如 (28) 所 示 的 问题 . 同样 地 , 根据 强 对 偶 性 , 式 (28) 可 
以 等 效 写成 (29). 而 (29) 同 样 可 以 由 式 (27) 表 示 . 因此 ， 式 (27) 对 于 任意 五 严格 正定 或 者 半 正 定 的 情 
况 都 是 准确 和 适用 的 . 


min (h” +v; V)E + ho — Q; (w,) >0 (28) 
ES 
dualizing 
max hoi (ui) > 0. (29) 


hj (u;)=0,u1;,u2i 20 
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利用 半 正 定 约束 (27)， 问 题 (24) 可 以 精确 地 转化 成 一 个 半 正 定 规划 : 


A) ey a PEO) eh acre (30) 
s.t. (27) WieT. 
可 以 看 到 问题 (30) 具 有 有 限 多 个 以 问题 (23) 的 可 行 域 顶 点 为 指标 的 约束 . 为 了 处 理 这 个 问题 ， 我 们 采 
用 延迟 约束 生成 法 来 求解 问题 (30) ， 具 体 过 程 如 下 所 述 . 
首先 ， 给 定 顶 点 的 子 集 (OAc pz = yD}, IER 了 < 加， 求解 问题 (30) 的 松弛 问题 以 得 
到 一 个 待 校 验 的 解 (H4, h*, hi). 然后 ， 通 过 求解 如 下 的 双 目 问题 来 校 验 约束 (22b) 是 否 满足 [Wei et 


al.(2016)], 

R(H*,h*,hs) = min €'H*E +h E+ 

EAH WY 
(g'A— (Ax dj 1 一 (VE Ww) 'v) 
s.t. co-(G'A+B'w+U'v)=0 
ALO p>O,EEE. 

如 果 R(H*, h*, ho) 之 0, 则 约束 (22b) 是 满足 的 ,可 以 确定 当前 解 是 问题 (22) 的 最 优 解 . 如 果 R(H*,h*, ho) > 
0 不 成 立 , 则 存在 一 组 对 偶 变量 的 解 A et} 将 其 记 为 {Ar Ar ri ， 这 个 解 是 一 个 新 的 
顶点 ， 可 以 将 其 扩充 到 问题 (30) 中 ， 再 次 求解 松弛 问题 并 重复 上 述 过 程 . 


(31) 


3.2 概率 分 布 生成 法 


用 上 文 得 到 的 半 正 定 规划 (30) 替 换 (12) 的 第 二 阶段 问题 ， 可 以 得 到 如 下 的 混合 整数 半 正 定 规划 ， 


i 
min cr © + Z(a). (32) 


混合 整数 半 正 定 规划 问题 (32) 可 以 采用 延迟 约束 生成 法 的 框架 ， 通 过 序 贯 求解 多 个 大 规模 的 混合 整 
数 半 正定 规划 问题 来 求解 受 限于 现 有 混合 整数 羊 正定 规划 求解 器 的 能 力 ， 延 迟 约束 生成 法 实际 上 
是 不 适用 的 [Zheng et al.(2020)]. 

为 了 高 效 求解 分 布 鲁 棒 机 组 组 合 问题 , 可 以 借鉴 文献 [Zeng and Zhao(2013)] 的 作者 提出 的 原 割 
平面 算法 ( 即 列 和 约束 生成 法 )， 该 算法 可 以 非常 有 效 地 解决 两 阶段 鲁 棒 优化 问题 . 列 和 约束 生成 法 
的 基本 思想 是 对 于 迭代 过 程 中 每 一 个 辨识 到 的 极端 场景 ,在 原 变量 的 决策 空间 不 断 地 构造 运行 约束 ， 
使 得 第 一 阶段 问题 满足 极端 场景 的 运行 要 求 . 本 文 在 列 和 约束 生成 法 的 框架 下 ， 提 出 一 种 原 割 平面 
算法 来 有 效 处 理 分 布 鲁 棒 机 组 组 合 问题 . 考虑 到 该 算法 是 基于 迭代 过 程 中 辨识 到 的 极端 概率 分 布 而 
非 极端 场景 ， 后 续 我 们 称 其 为 概率 分 布 生成 法 . 

下 文 我 们 首先 揭示 一 个 重要 的 原理 和 结论 ， 即 求解 半 正定 规划 问题 (30) 之 后 可 以 得 到 特定 的 第 
一 阶段 决策 z 相应 的 极端 概率 分 布 (极限 分 布 )， 然 后 再 介绍 求解 所 提 的 分 布 鲁 棒 机 组 组 合 问题 求 
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对 问题 (30) 再 次 使 用 对 偶 方 法 ， 可 以 得 到 原 问 题 (21) 有 限 维 形式 的 版 本 : 


Z(x)= max 》 piQ; (x, £i) — (Wm) vi (33a) 
Oi, Ti,Pi i 
si S p= (33b) 
icT 
yo mae (33c) 
i€I 
Soo, 3 B+ EE" (33d) 
i€I 
Oi Ti 
~>0 VicT (33e) 
T; pi 
pE™ < m; <p ViET, (33f) 


其 中 变量 oj, mi, pi 是 如 算法 1 的 第 1 行 所 示 的 对 偶 变 量 . 不 失 一 般 性 ， 这 里 假设 p; 的 取 值 均 为 非 
零 值 ， 从 约束 (33b)-(33c) 可 以 看 出 ， 变 量 的 组 合 表达 式 zi/p; (后 文 记 为 é) 是 与 随机 向 量 紧 密 相 
FRAY. 具体 地 ， 第 i 个 标量 pi 表示 事件 &; 的 概率 质量 .约束 (33e) 等 效 表 征 两 个 约束 ， 即 pi > 0 
Al mim! /pi = pi&i&i < oi. 因此， 约束 (33b)-(33e) 表 明 问 题 (30) 的 对 偶 变 量 构成 了 一 个 离散 分 布 
{(é&i, pi)li 站， 该 分 布 的 期 望 和 协 方差 分 别 等 于 和 小 于 E. 进一步 地 ， 约 束 (33f) 还 表明 该 离散 分 
布 的 支撑 集 是 集合 E. 根据 上 述 分 析 可 知 ， 离 散 分 布 属 于 模糊 集 (14). 基于 强 对 偶 性 ， 对 偶 变 量 的 最 
优 值 使 得 问题 (33) 取 得 最 大 值 以 及 问题 (30) 的 最 小 值 ， 所 以 所 构建 的 概率 分 布 是 使 问题 (21) 在 给 定 第 
一 阶段 方案 x 之 下 取得 最 大 值 的 极限 分 布 . 


半 正 定 规划 问题 隐 含 的 极限 分 布 可 以 在 搜索 第 一 阶段 问题 最 优 解 方 面 发 挥 关键 作用 . 与 标准 的 
列 和 约束 生成 法 类 似 ， 可 将 每 一 步 欠 代 辩 识 到 的 极限 分 布 扩充 到 第 一 阶段 问题 ， 对 第 一 阶段 问题 的 
可 行 域 形成 附加 的 约束 ， 从 而 得 到 更 接近 原 问 题 要 求 而 的 解 ， 循 环 直至 得 到 原 问题 的 最 优 解 . 


总 体 的 求解 方法 如 算法 1 所 示 . 第 ] 行 到 第 ] 行 是 主 循环 ， 该 部 分 迭代 求解 混合 整数 线性 规划 主 问 
题 和 半 正 定 规划 子 问题 ， 第 1 行 所 示 的 主 问题 采用 极限 分 布 和 Q, E) 的 原 问题 形式 来 生成 原 割 平 
面 . 每 个 主 循环 中 ， 需 要 求解 多 个 半 正 定 规划 子 问题 ， 直 到 约束 (22b) 满 足 ， 即 R(H*,h*, hg) > 0. 
需要 注意 半 正 定 规划 中 的 hos (ui) 采用 了 当前 第 一 阶段 问题 的 最 优 解 zxx. 在 第 1 行 ， 调 用 了 交替 优 
化 算法 来 求解 双 凸 问题 (31), 该 算法 轮流 地 固定 {A, 4,v} 来 求解 二 次 规划 问题 ,以 及 固定 & 来 求解 
线性 规划 问题 [Wei et al.(2016), Zheng et al.(2020)]. 交 蔡 优化 算法 是 求解 双 凸 规划 问题 的 局 部 最 优 
算法 ， 本 文采 用 多 个 初始 解 来 保证 算法 所 得 解 的 质量 .我们 仔细 校 验 了 交替 优化 算法 的 性 能 ， 该 算 
法 求解 本 文 的 双 凸 规划 问题 得 到 的 解 的 最 优 性 与 内 点 法 求解 器 IPOPT 相 比 是 具有 竞争 力 的 ， 而 在 
计算 效率 方面 则 显著 优 于 IPOPT. 尽管 算法 1 包含 了 启发 式 过程 ( 即 交替 优化 算法 )， 在 算 例 测试 中 
算法 1 的 收敛 性 仍然 是 非常 优良 的 ， 且 其 最 终 可 以 得 到 高 质量 的 开 停机 方案 . 
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算法 1 分 布 鲁 棒 机 组 组 合 问题 的 总 体 求解 方法 
o 选择 收敛 容 差 o 记 变量 r 的 最 优 值 为 as z*; 4 UB=00, LB=-%; 令 T=1, K =0; 
& {A1, 1, v1} 为 初始 顶点 . 
0: while UB — LB > e do 
利用 已 有 的 割 平 面 ， 求 解 如 下 主 问 题 以 获得 08, {MILP} 


v= min cilz+0 
LEX A>—00 
Ik 
s.t. 0 > XO pei Q(x, mii/ Pri) k= 1, wo, K. 
i=1 


o KeK+1, axe at, DBOvVUDIg cl. 
0: repeat 


0: 给 定 zx， 求解 如 下 子 问 题 以 获得 9, {SDP} 
9=crzrt, min tr(H'O) +h E+ ho 


=0,h,ho 1: Uzi 


H d =0: k "| 


1/2h;(w;)" hoi (ui) l m; Pi 


Uii, U2; > 0 t= 1, «te: 


s.t. 


0: if R(H*,h*,hġ) < 0 then {QP 和 LP 迭代 } 
0: Ik 4} Ik +1, {Ars Mie Vig} 全 TAN Dr 
0: end if 
0: until R(H*,h*,hġ) > 0 
0: (Tki PK i) < (T7, PF) i=1,..IK, UB 1. 
0: end while 
0: 返回 TK. =0 
4 算 例 研究 


本 节 开 展 数 值 实验 来 验证 风 功 率 预 测 的 条 件 误差 估计 器 的 有 效 性 ， 以 及 证 明 所 提出 的 分 布 鲁 棒 
机 组 组 合 模型 及 其 求解 方法 的 有 效 性 . 实验 采用 了 4 个 修改 的 IEEE 标准 测试 系统 . 如 表 1 所 示 ， 测 
试 系统 的 节点 数量 小 至 6 个 节点 ， 大 至 300 个 节点 ， 其 中 300 节点 系统 有 8 个 风电 场 ， 每 个 系统 
的 最 大 风 功 率 占 比 均 达 到 15%. 本 文采 用 的 机 组 组 合 模型 ， 以 及 测试 系统 的 完整 数据 均 可 以 在 线 获 
取 [Zheng(2020)]. 

实验 中 的 MILPs、LPs 和 QPs 均 采 用 GUROBI 8.1 求解 ，SDPs 采用 MOSEK 9.1 求解 . 求 
解 器 经 YALMIP 调用 . 除非 特殊 说 明 ， 混 合 整数 的 相对 收敛 间隙 默认 设置 为 10“， 连 续 优 化 问题 的 
收敛 容 差 都 采用 YALMIP 的 默认 值 . 算法 1 的 。 值 设 置 为 . 所 有 计算 均 在 一 台 运 行 频率 为 3.00GHz、 
内 存 为 8GB、CPU X Intel i5-8250U 的 个 人 计算 机 上 完成 . 


4.1 基于 Lasso 的 估计 器 的 有 效 性 


实验 的 风 功 率 数据 来 源 于 泛 加 拿 大 风电 并 网 研究 的 网 站 开源 数据 ,https://canwea.ca/wind-integration-study/ 
wind-data/. 实验 利用 了 3 年 的 日 前 预测 和 实时 数据 ， 整 个 数据 集 的 前 一 半 用 于 统计 ， 后 一 半 用 于 
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表 1 测试 系统 的 配置 


系统 6 节点 系统 ”39 节点 系统 118 节点 系统 300 节点 系统 
发 电机 数量 4 10 54 69 
风电 场 数量 1 2 3 8 
总 装机 容 = (MW) 630 7367 7220 32678 
jel ee ae 量 (MW) 110 1300 1274 5767 
150 
-=Lasso 
一 无 偏 估计 
入 140 -ARH 
= 
= 
~ 130 
jay 
(ep) 
三 
æ 120 


1 5 9 13 17 2 25 29 33 37 41 45 48 
风电 变量 序号 


图 1 风 功 率 预测 误差 均 方 根 的 对 比 . 


仿真 测试 . 

我 们 将 本 文 提出 的 基于 Lasso 的 条 件 误 差 估计 器 的 实际 效果 与 无 偏 窍 估计 器 的 效果 进行 对 比 ， 
后 者 对 风 功 率 预测 误差 的 估计 值 是 历史 误差 的 平均 值 ， 即 € = 48 E. 图 1 给 出 了 各 项 预测 结 
的 均 方 根 (RMS) 指标 ， 包 括 日 前 预测 值 ( 记 为 “ 原 曲线 ”)、 经 由 基于 Lasso 的 条 件 误差 估计 校正 
后 的 预测 曲线 〈 记 为 “Lasso”) 、 经 由 无 偏 矩 估计 校正 后 的 预测 曲线 〈 记 为 “无 偏 估计 ”) . 

在 这 个 案例 中 ,39 节点 系统 包含 了 两 个 650MW 的 风电 场 ,因此 随机 向 量 & 的 维度 是 48. 将 条 件 
误差 的 估计 值 苹 加 到 日 前 预测 值 后 ， 实 际 预测 误差 的 均 方 根 有 所 降低 . 均 方 根 最 大 降低 了 15.94MW 
( 即 功 率 序列 上 的 第 25 个 值 ) . 如 图 1 所 示 ， 无 偏 估计 器 的 误差 估计 值 准确 的 不 高 ， 对 日 前 预测 起 不 
到 精度 提升 的 作用 . 

由 于 在 计算 日 前 机 组 组 合 问题 的 时 候 ， 前 一 天 的 日 前 预测 值 和 实际 值 大 部 分 是 已 知 的 ， 可 以 将 
前 一 天 的 日 前 预测 值 和 实际 值 作为 Lasso 的 附加 输入 . 将 考虑 了 附加 输入 的 模型 记 为 Lasso*. 实际 
o Lasso* 之 后 RMS 会 进一步 降低 . 然而 ， 当 系统 存在 多 个 风电 场 时 ，Lasso” 模型 的 训 


练 会 变 得 比较 耗 时 . 例如 ， 


在 含有 8 个 风电 场 的 案例 中 ， 训 练 数据 集 Å 的 规模 为 548 x 576 (数据 


集 含 548 个 样本 ， 每 个 样本 的 维 


度 是 24 x 8 x 3) 


. 如 果 在 100 个 候选 的 正则 化 参数 入 中 通过 2 层 


的 交叉 检验 来 选择 最 


佳 参 数 ， 并 将 坐标 下 降 算法 的 迭代 数 上 限 设置 为 10”， 


则 在 192 个 Lasso 模型 


需要 耗费 超过 2 小 时 来 训练 每 个 模型 . 为 了 节省 训练 时 间 ， 对 于 含 8 个 风电 场 的 案例 ， 我 们 将 和 参 
数 候选 值 的 数量 减少 到 10 个 ， 将 和 迭代 数 的 上 限 降低 为 104. 通过 这 一 措施 ， 该 案例 的 所 有 Lasso 模 
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表 2 不 同方 法 应 用 在 不 同 风电 场 集合 上 所 取得 的 风 功 率 均 方 根 误差 
风电 场 数量 Lasso% Lasso ”无 偏 估计 原 曲 线 
RMS (MW) | 17.79 18.93 19.31 19.25 
改进 程度 7.56% 1.63% -0.32% - 
RMS (MW) | 17.91 18.78 19.57 19.49 
改进 程度 8.10% 3.62% -0.41% - 
RMS (MW) 18.02 18.90 19.67 19.59 


3 
改进 程度 8.03% 3.54% -0.40% > 
RMS (MW) 17.35 17.90 18.29 18.44 


改进 5.89% 2.93% 0.81% = 


型 的 训练 时 间 降 低 到 了 10 个 小 时 ， 即 每 个 模型 耗费 大 约 3 分 钟 . 
> 表 2 给 出 了 经 由 Lasso*、Lasso、 无 偏 窍 估计 器 等 模型 校正 后 和 原 预测 值 相对 应 的 风 功 率 预 测 误 
Al 差 的 均 方 根 . 所 有 风电 场 的 容量 均 归 算 成 了 LOOMW. 可 以 看 出 原 预 测 值 的 均 方 根 误差 约 为 20%. 使 
=) 用 基于 Lasso 的 条 件 误差 估计 器 进行 校正 之 后 ， 风 功率 的 预测 精度 最 大 可 以 提升 8% (Lasso*) ， 进 
而 可 以 进一步 得 到 更 准确 的 模糊 集 ， 从 表 2 还 可 以 看 出 ， 在 含 8 个 风电 场 的 案例 中 Lasso 降低 均 方 
根 误 差 的 幅度 较 小 ， 即 为 了 提升 训练 效率 会 牺牲 校正 精度 . 实际 上 ，Lasso 模型 可 以 离线 进行 并 行 训 
练 ， 所 以 可 以 采用 更 精细 的 训练 方式 以 获得 更 高 的 校正 精度 . 
Č 我 们 用 相对 Frobenius 距离 来 评估 协 方差 矩阵 的 估计 精度 . 将 估计 的 协 方差 矩阵 记 为 S, 从 样本 
© 外 数据 计算 得 到 的 协 方差 矩阵 记 为 >o, 则 相对 Frobenius 距离 的 表达 式 为 Vtr((D — Xo)( — Xo)! )/ tr (P00 ). 
N 在 4 个 案例 中 ， 相 对 Frobenius 距离 的 值 分 别 是 29.52%, 28.81%, 29.52% 和 31.66%. 


4.2 ” 开 停 机 计划 的 经 济 性 和 和 鲁 棒 性 


= WEIR, 7 2]— BAB hitter Za, Xa ACSI BW. 为 了 更 准确 模 
拟 两 阶段 的 运行 模式 (即日 前 机 组 组 合 、 日 内 经 济 调度 ) ， 逐 日 进行 仿真 . 对 于 数据 集中 的 每 一 天 ， 
基于 日 前 功率 预测 计算 预测 误差 的 条 件 期 望 ， 得 到 的 条 件 期 望 和 协 方差 矩阵 用 于 构造 模糊 集 ， 基于 
模糊 集 可 以 求解 分 布 鲁 棒 机 组 组 合 问 题 ， 获 取 开 停机 方案 . 最 后 ， 利 用 当天 的 实际 风 功 率 数据 ， 计 
算 给 定 开 停机 方案 所 对 应 的 真实 运行 费用 . 

本 文采 用 了 线性 的 费用 函数 . 弃 风 和 负荷 中 断 的 价格 分 别 设置 为 2008/MW 和 1,0008/MW [Wei 
et al.(2016)]. 为 了 计 及 系统 调节 能 力 不 足 带 来 的 高 额 惩罚 费用 ， 将 模型 的 松弛 变量 的 惩罚 系数 设置 
为 5,000$/MW [Bertsimas et al.(2013)]. 在 实际 实验 结果 中 ， 松 弛 变量 的 优化 值 均 为 0. 

我 们 将 从 本 文 所 提 的 分 布 鲁 棒 优 化 模型 (后 续 记 为 模型 1) 得 到 的 运行 费用 与 其 他 三 个 模型 的 
费用 进行 对 比 , 包括 : 不 考虑 风 功 率 预 测 的 条 件 误 差 和 误差 相关 性 的 模型 2， 即 模糊 集 (14) 的 第 二 行 
WE ERO, 第 三 行 不 考虑 ; 基于 数据 驱动 的 随机 规划 的 模型 3, 该 模型 使 用 [Dupacova et al.(2003)] 
提出 的 基于 距离 的 场景 缩减 算法 来 从 训练 数据 集 (FF 548 个 场景 ) 生成 50 个 场景 9?; 文献 [Veloso 
et al.(2019)] 提出 的 基于 数据 驱动 的 鲁 棒 优化 的 , 该 模型 的 不 确定 集 由 历史 发 电 时 序数 据 集 的 凸 包 构 
成 , 因此 不 需要 整定 不 确定 集 参 数 (如 不 确定 性 预算 [Bertsimas et al.(2013), Wei et al.(2016), Zheng 


1) 缩减 算法 可 在 线 获取 : https://www.mathworks.com/matlabcentral/fileexchange/67909-scenred. 
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本 ene PSE Sm 
T 一 一 一 一 ”| 5 
548 天 的 日 前 预测 和 
实时 数据 估计 风 功 率 ESE EUS 
六 一 5 预测 到 实时 误差 的 双 
的 矩 Lasso “上 H ENSA 
1， 条 人 | mpiya 
$; =s $ 望 
二 1096 天 的 | | HAD - 


anne 利用 风 功率 的 实时 数据 逐日 运行 
开 停机 方案 | 进行 实时 经 济 调度 模拟 ”| 运行 费用 


图 2 数值 实验 的 框架 . 


et al.(2020)]) . 这 四 个 模型 均 为 数据 驱动 的 模型 4， 可 以 使 得 对 比 更 加 公平 此外， 模型 运行 费用 
的 对 比 采取 了 图 2 所 示 的 逐日 仿真 框 染 . 

为 了 节省 仿真 时 间 ，39 节点 、118 节点 、300 节点 系统 只 考虑 了 关键 输电 线路 的 传输 能 力 约束 . 
尽管 目前 学 术 界 已 经 有 一 些 减轻 输电 线路 安全 约束 带 来 的 计算 负担 的 精细 方法 ， 但 不 失 一 般 性 ， 我 
们 采取 了 减少 建 模 线路 的 数量 的 简单 方法 . 具体 地 ， 首 先 求解 一 个 确定 性 的 安全 约束 机 组 组 合 问题 
实例 , 然后 从 中 挑选 传输 功率 最 接近 输电 容量 的 10 条 线路 ,后续 仅 对 这 些 线路 建 模 . 此 外 , 从 第 549 
天 和 第 1096 天 的 所 有 样本 外 数据 中 选 出 20% 的 天 数 进行 仿真 . 感 兴趣 的 读者 可 以 从 [Zheng(2020)] 
了 人 解 测试 系统 配置 和 仿真 实验 的 更 多 细节 . 

表 3 总 结 了 运行 费用 .对 于 每 个 测试 系统 : 定义 了 平均 费用 ， 表 示 所 有 仿真 日 的 平均 运行 费用 ; 
定义 了 费用 差异 ， 表 示 该 模型 的 每 天 的 费用 与 4 个 模型 当天 平均 费用 的 差异 值 . 因此 ， 平 均 差 异 反 
映 了 每 个 模型 的 费用 节省 程度 、 最 小 差异 和 最 大 差异 分 别 反映 了 每 个 模型 的 最 大 费用 节省 量 和 最 大 
费用 多 消耗 量 . 

从 表 3 可 以 看 出 ， 除 了 在 6 节点 系统 外 (在 这 个 系统 是 模型 4 具有 最 小 的 平均 费用 ) 模型 1 的 


平均 费用 是 最 低 的 . 在 118 节点 系统 ， 模 型 1 的 平均 费用 节省 量 相 比 模型 4 的 高 达 1.20%. 模型 1 
的 最 大 费用 节省 量 总 是 最 大 的 ， 可 以 高 达 平 均 费 用 的 (如 模型 1 在 118 节点 系统 的 “最 小 差异 ” 行 


Bras); 模型 1 的 “最 大 差异 ”项 也 比较 小 ， 这 意味 着 其 开 停机 方案 在 费用 方面 的 鲁 棒 性 更 强 . 总 体 
而 言 ， 基 于 数据 驱动 的 随机 规划 的 模型 3 的 性 能 是 最 为 平均 的 ， 而 基于 数据 驱动 的 鲁 棱 优化 的 模型 
4 的 性 能 波动 较 大 . 所 以 ， 按 照 模 型 1 的 思路 基于 极限 分 布 制定 的 开 停 机 方案 ， 可 以 比 模型 3 和 模 
型 4 更 好 地 平衡 运行 的 经 济 性 和 重 棒 性 . 


4.3 ”极限 分 布 的 特征 及 其 对 开 停 机 计划 的 影响 


在 模型 工 的 模糊 集 利 用 风 功 率 预 测 误差 的 条 件 期 望 提 升 了 预测 精度 ， 利 用 协 方差 矩阵 来 排除 不 
实际 的 概率 分 布 . 对 比 模型 1 和 模型 2 可 以 发 现 ， 分 布 鲁 棒 模型 得 到 的 优化 方案 的 经 济 性 取决 于 模 
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表 3 运行 费用 的 对 比 (k$) 
系统 费用 模型 1 模型 2 模型 3 模型 4 
平均 费用 | 162.87 162.96 162.91 162.83 
Se 平均 差异 | -0.02 0.07 0.02 -0.06 
6 节点 系统 | a a 
最 小 差异 | -1.20 -0.76 -1.00 -0.42 
最 大 差异 | 0.19 0.45 0.17 2.97 
平均 费用 | 1,254.28 1,254.72 1,266.16 1,263.83 
ibe 平均 差异 | -5.47 -5.03 6.41 4.08 
39 节点 系统 | -人 
最 小 差异 | -53.17 -53.17 -48.62 -5.03 
最 大 差异 | 0.42 0.86 14.56 154.96 
平均 费用 | 1,608.54 1,609.21 1,616.20 1,628.11 
平均 差异 | -6.97 -6.30 0.68 12.59 
118 节点 系统 平均 差异 | -6.97 
最 小 差异 | -84.34 -83.52 -77.38 -8.30 
最 大 差异 | 4.19 3.92 9.48 245.25 
平均 费用 | 3,484.28 3,500.02 3,496.85 3,495.55 
平均 差异 | -9.90 5.84 2.68 1.38 
300 节点 系统 | PPT) 9.90 
最 小 差异 | -37.55 -8.43 -8.53 -18.51 
最 大 差异 | 0.61 23.34 16.75 24.94 


糊 集 或 者 极限 分 布 的 质量 . 


图 3 


反 虑 条 件 误差 和 误差 相关 性 得 到 的 极 


限 分 布 . 


如 第 3.2 节 所 述 ， 子 问题 对 偶 变 量 的 最 优 值 可 以 构造 极限 分 布 {(&i, pi)li E Zh, RP &; 代表 第 i 
个 事件 ，p; 代表 该 事件 的 概率 . 显然 ， 模 型 1 和 模型 2 得 到 的 极限 分 布 是 不 一 样 的 . 我 们 有 必要 将 
这 些 极限 分 布 进 行 可 视 化 ， 以 便 直观 地 分 析 极 限 分 布 的 形态 是 
出 了 39 节点 系统 某 个 实例 中 的 极限 分 布 . 
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图 4 不 考虑 条 件 误差 和 误差 相关 性 得 到 的 极限 分 布 . 


为 了 使 展示 的 结果 更 加 人 简洁 清晰 ， 每 个 图 以 4 个 子 图 的 形式 给 出 分 布 中 的 概率 较 高 的 前 8 个 事 
件 . 图 中 实 线 代表 分 布 中 的 事件 ， 线 的 粗细 则 表示 了 该 事件 概率 的 大 小 ， 图 3 和 图 4 用 带 框 标记 的 线 
画 出 了 【校正 后 的 ) 日 前 风 功 率 预测 曲线 . 可 以 看 出 模型 1 的 极限 分 布 排除 了 很 多 频繁 和 快速 假 坡 
的 风 功 率 场景 ， 因 为 这 些 场景 会 导致 较 大 的 协 方差 . 尽管 图 4 的 事件 是 以 期 望 值 为 中 心 ， 由 于 没有 考 
虚 时 空 关联 性 ， 其 事件 的 时 空 波动 过 于 剧烈 ,值得 注意 的 是 ， 方 差 参 数 是 不 足以 用 于 对 时 空 波 动 和 
关联 性 进行 建 模 的 [Babaei et al.(2019)]， 因 此 有 必要 在 构建 模糊 集 的 时 候 考虑 完整 的 协 方差 矩阵 . 

为 了 应 对 风 功 率 的 快速 疏 坡 事件 ， 需 要 多 开 灵 活 和 昂贵 的 发 电机 组 ， 且 机 组 的 计划 会 出 现 更 加 
频繁 的 开 停 机 切换 . 所以， 排除 这 些 不 实际 的 极端 事件 对 于 降低 开 停 机 计划 的 保守 性 是 尤其 关键 的 . 
在 如 上 所 示 的 39 节点 系统 中 ， 模 型 1 和 模型 2 所 得 的 10 台 机 组 的 开 停机 计划 如 图 5 所 示 . 可 以 看 
出 机 组 4 (变量 序号 从 73 到 96) 和 机 组 5 (变量 序号 从 97 到 120) 基于 模型 2 是 被 调用 的 ， 且 机 
组 5 的 启 停 更 加 频繁 . 尽管 两 个 模型 的 开 停机 方案 均 不 会 引发 实际 运行 中 的 切 负荷 和 弃 风 ， 模 型 2 
的 运行 费用 是 显著 高 于 模型 1 的 (主要 是 因为 额外 的 空 载 成 本 和 启 停 成 本 ) . 所 以 ,模型 2 的 开 停 
机 方案 是 过 于 保守 的 . 这 进一步 印证 了 条 件 误差 和 协 方差 信息 对 分 布 鲁 棒 机 组 组 合 模型 决策 质量 的 
重要 影响 . 


4.4 ”计算 效率 


表 4 给 出 了 求解 每 个 模型 的 耗 时 与 闪 代 次 数 . 对 于 118 节点 和 300 节点 系统 ， 主 问题 的 混合 整 
数 收银 间隙 设置 为 10°. 对 于 300 节点 系统 ， 如 果 运 行 时 间 超过 2 个 小 时 则 耗 时 数据 是 空 的 ; 对 于 
其 他 系统 ， 如 果 运 行 时 间 超 过 1 个 小 时 则 耗 时 数据 是 空 的 . 考虑 到 Lasso 模型 的 训练 是 可 以 离线 进 
行 的 ， 且 不 需要 每 天 更 新 训练 ， 耗 时 项 并 不 包括 训练 时 间 . 

对 于 除 模型 3 之 外 (模型 3 是 一 个 MILP 问题 ， 可 用 求解 器 一 次 性 求解 ) 的 所 有 模型 ， 均 记录 
了 求解 所 需要 的 迭代 次 数 ， 并 在 耗 时 数据 的 后 方 列 出 . 模型 4 是 容易 求解 的 . 使 用 列 和 约束 生成 法 
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开 停 机 变量 的 序号 


5 39 节点 系统 案例 的 开 停机 方案 对 比 . 


求解 数据 驱动 的 鲁 棒 优 化 模型 ， 所 需要 的 迭代 次 数 均 在 T 次 以 内 ， 耗 时 均 在 3 分 钟 之 内 . 因为 考虑 
了 50 个 场景 ， 问 题 规模 较 大 ， 求 解 模型 3 的 计算 负担 稍微 高 于 模型 4. 模型 1 和 模型 2 的 计算 负 
担 明显 更 大 . 在 大 系统 的 测试 案例 中 ， 模 型 1 比 其 他 模型 需要 更 多 的 计算 时 间 . 这 有 两 方面 的 原因 : 
i) 半 正 定 规划 子 问题 的 规模 随 着 随机 向 量 的 维 数 升 高 而 升 高 ， 而 目前 半 正 定 规划 求解 需 的 能 力 仍 比 
较 有 限 ; ii) 混合 整数 线性 规划 主 问题 的 规模 随 着 极限 分 布 不 断 添加 而 增 大 . 


在 300 节点 系统 的 实例 中 ， 需 要 求解 多 个 大 规模 的 半 正 定 规划 问题 : 由 于 系统 有 8 个 风电 场 ， 
随机 向 量 的 维度 是 192 x 1， 半 正定 规划 的 决策 变量 由 27 个 192 x 1 维 向 量 和 了 个 193 x 193 AER 
阵 组 成 ， 其 中 工 是 当前 子 问题 的 序号 . 虽然 理论 上 半 正 定 规划 可 在 多 项 式 时 间 内 求解 ， 但 是 当前 基 
于 内 点 法 的 求解 器 仍然 无 法 有 效 处 理 超 大 规模 的 半 正 定 规划 问题 . 因此 ， 为 了 在 合理 与 可 控 的 时 间 
内 (如 日 前 现货 市 场 出 清 的 4 小 时 时 间 限 制 ) ， 我 们 将 300 节点 系统 的 测试 案例 的 子 问 题 项 点 数 限 
制 在 5 以 内 . 通过 这 个 设置 ， 子 问题 将 返回 一 个 事件 数量 不 超过 5 的 可 行 概率 分 布 (不 一 定 是 极限 
分 布 ). 实际 计算 结果 表明 ， 返 回 的 可 行 概率 分 布 仍然 是 比较 苛刻 的 分 布 ， 且 算法 1 能 够 收敛 到 高 质 
量 的 开 停机 方案 ( 详 见 第 4.2 节 的 表 3) ， 


这 一 节 也 将 算法 1 与 其 他 两 个 求解 方法 相对 比 , 即 延 迟 约束 生成 法 和 Benders 方法 . 文献 [Zheng 
et al.(2020)] 采用 延迟 约束 生成 法 , 在 儿 分 钟 内 成 功 求解 了 6 节点 系统 的 分 布 鲁 棒 机 组 组 合 问题 ， 然 
而 该 模型 的 支撑 集 是 无 界 的 ， 与 本 文 的 模型 不 同 ， 如 果 使 用 延迟 约束 生成 法 将 问题 (32) 作 为 混合 整 
数 半 正 定 规划 直接 求解 ,或 者 使 用 Benders 方法 ， 则 在 6 节点 系统 的 案例 中 算法 也 无 法 于 1 小 时 内 
收敛 , 表格 也 给 出 了 这 两 个 求解 方法 在 求解 时 间 达 到 上 限时 典型 的 迭代 次 数 . 根据 表 4， 算 法 1 在 求 
解 两 阶段 分 布 鲁 棒 机 组 组 合 问题 方面 的 能 力 是 明显 优 于 其 他 两 个 对 比方 法 的 . 
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表 4 不 同 模型 或 者 方法 的 耗 时 和 迭代 次 数秒 、 次 ) 


系统 模型 1 
Algorithm 1 MISDP Benders 
6 节点 系统 | 80.60 (1-5) - (41) - (40) 
39 节点 系统 | 217.24 (2-5) - (17) s7) 
118 节点 系统 | 2348.21 (2-4) - (7) - (3) 
300 节点 系统 | 5142.60 (2-3) - (2) - (2) 
系统 模型 2 模型 3 模型 4 


6 节点 系统 13.33 (1 
39 节点 系统 | 44.29 (2- 
118 节点 系统 | 1091.62 
300 节点 系统 | 2111.53 


( 
) 163.43 59.34 (2- 
) 714.84 189.33 (2-7) 


5 结论 与 讨论 


基于 Lasso 的 条 件 误差 估计 器 有 效 提升 了 风 功 率 预 测 的 精度 . 对 风 功 率 条 件 预测 误差 和 误差 相 
关 性 信息 的 利用 ， 有 效 地 排除 了 不 符合 实际 情况 的 风 功 率 概率 分 布 ， 因 此 得 到 的 分 布 鲁 棒 机 组 组 合 
模型 能 够 确定 更 经 济 和 可 靠 的 方案 ， 本文 给 出 的 混合 整数 半 正 定 规划 模型 是 含 二 阶 矩 约束 和 多 面体 
支撑 集 的 分 布 鲁 棒 机 组 组 合 模型 的 等 价 模型 ， 所 提出 的 极限 分 布 生成 方法 相 比 已 有 方法 能 够 更 加 高 
效 地 求解 两 阶段 分 布 鲁 棒 机 组 组 合 问题 . 下 一 步 工 作 将 致力 于 提升 二 阶 矩 估计 的 精度 ， 提 出 其 他 类 
型 两 阶段 分 布 鲁 棒 优化 问题 的 求解 算法 ， 以 及 利用 半 正 定 规划 的 特殊 结构 提升 算法 效率 等 方面 . 
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Abstract In this paper, a study of the day-ahead unit commitment problem with stochastic wind power genera- 
tion is presented, which considers conditional and correlated wind power forecast errors through a distributionally 
robust optimization approach. Firstly, to capture the characteristics of random wind power forecast errors, the 
least absolute shrinkage and selection operator (Lasso) is utilized to develop a robust conditional error estimator, 
while an unbiased estimator is used to obtain the covariance matrix. The conditional error and the covariance 
matrix are then used to construct an enhanced ambiguity set. Secondly, we develop an equivalent mixed integer 
semidefinite programming (MISDP) formulation of the two-stage distributionally robust unit commitment model 
with a polyhedral support of random variables. Further, to efficiently solve this problem, a novel cutting plane 
algorithm that makes use of the extremal distributions identified from the second-stage semidefinite programming 
(SDP) problems is introduced. Finally, numerical case studies show the advantage of the proposed model in cap- 
turing the spatiotemporal correlation in wind power generation, as well as the economic efficiency and robustness 
of dispatch decisions. 


Keywords Column-and-constraint generation, distributionally robust optimization, extremal distribution, spa- 


tiotemporal correlation, wind power forecast. 


21 


